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It is demonstrated theoretically that the counter-rotating quadrupole mode in a vortex of Bose- 
Einstein condensates can decay into a pair of Kelvin modes via Beliaev process. We calculate the 
spectral weight of density-response function within Bogoliubov framework, taking account o f both 



Beliaev and L andau processes. Good agreement with experiment on Rb by Bretin et al. [cond 



mat/0211101] allows us to unambigiously identify the decayed mode as the Kelvin wave propagating 
along a vortex line. 
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Much attention has been focused on Bose-Einstein con- 
densation (BEC) realized in alkali-atom gases fl. Quan- 
tized vortices which are a hallmark of superfluidity 
are created in various methods, such as phase imprinting 
||, mechanical rotation by optical spoon and topo- 
logical Berry phase engineering ^ . Hundreds of vortices 
are trapped in a BEC system Q . The creation and decay 
processes of vortices are investigated experimentally and 
theoretically, giving rise to a general consensus Q that 
quantized vortices in a scalar, i.e. one-component BEC 
are well described by the Bogoliubov framework regard- 
ing to the static properties, such as the density profile or 
core radius, etc. 

In contrast, regarding to the dynamical aspects, the 
study of low- lying collective modes is rather scarce in the- 
ory and particularly in experiment. Needless to say, the 
low-lying Fermionic excitations in a vortex have played a 
fundamental role in charged or neutral Fermion systems, 
that is, the mixed state in a superconductor and su- 
perfluid 3 He. Here we have a unique opportunity to in- 
vestigate Bosonic excitations associated with a vortex, 
which was difficult in superconductivity. In particular 
the so-called Kelvin mode || propagating along the vor- 
tex line, which is studied in classical normal fluids and 
superfluid 4 He, is interesting to identify and characterize 
in the present dilute Bose gases. This is an unexplored 
region. 

Note that Bosonic excitations with lower energy in a 
vortex-free BEC are thoroughly studied; the breathing or 
monopole mode with the azimuthal angular momentum 
qe = 0, the dipole Kohn mode qg = 1, and the quadrupole 
mode qe = 2 for an axis-symmetric system jl| . 

Recently, Bretin et al. || have performed an experi- 
ment to examine the quadrupole modes with qg = ±2 for 
a long-cigar shaped BEC with a vortex line, observing the 
decay process. Their results are summarized as follows: 
(1) The one of the splitted quadrupole mode qe = —2, 
which rotates opposite to the vortex winding, decays two- 
time faster than the other co-rotating quadrupole mode 
with qe = +2. (2) After the decay of qe = —2 mode, there 
remains a density oscillation pattern along the long axis 



(z-axis) whose nodes are 7 or 8 within the length of the 
condensate. The oscillation pattern is localized near the 
vortex core, seen in the radial direction profile. 

Here we investigate the physical implication of these in- 
teresting observations, by calculating the density-density 
response function based on the wave functions and eigen- 
values of the Bogoliubov-de Gennes equation for describ- 
ing the collective modes of Bosonic excitations. 

Before going into detailed computation of a cylindri- 
cal system with a vortex, we first give a clear physical 
picture for these phenomena which is fully justified mi- 
croscopically later: At low temperatures, among the two 
possible decay channels, the Beliaev process dominates 
over the remaining Landau process (lQ j. The counter- 
rotating quadrupole mode with the energy io-2 and the 
angular momentum qe = — 2 (in units of the radial har- 
monic frequency oj r and fi = 1 respectively) can de- 
cay into a pair of the dipole mode with conserv- 
ing the energy w_2 — > 2c<j_i and the angular momentum 
(qe = —2) — > (qg = —1) + (qe = —1). However, as seen 
from Fig. 1 where a Beliaev process is depicted as a func- 
tion of the wave number q z along the z-axis, the pair cre- 
ated dipole modes (see Fig. 1(b)) should have a finite q z , 
namely, (jj_2(g z = 0) ^u>-i(q z ) + uj-i(— q z ). Note that in 
the absence of a vortex shown in Fig. 1(a) there is no Be- 
liaev process because 2oji(q z ) > which is well-known, 
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FIG. 1: Solid circles show the dispersion relations u~\(q z ) of 
(qe, q r ) = (— 1, 0) along q z in the case of the vortex-free state 
(a) and a single- vortex state (b). ® denotes lu-2- 
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leaving only Landau decay active Likewise as for the 
other comoving quadrupole mode with qg = +2 there is 
no Beliaev process because 2w +1 > u 2 - Only the Landau 
process, in which the quadrupole mode annihilates with 
the thermally already excited modes into other modes, is 
responsible for the decay. 

Here it is central that the so-called anomalous mode 
uj-i at q z = has a negative values |l2|, [l3|, [l4| relative 
to the condensate at the zero energy and that their wave 
function is radially localized at the core region whose 
wave length is an order of the coherence length. Along 
the ^-direction it oscillates sinusoidally with q z . These 
explain the above experimental facts (1) and (2) simul- 
taneously. In addition the decay time of w+2 should be 
the same as in the vortex-free case, which is indeed the 
case ||. 

Let us now consider this novel phenomenon of 
the single-vortex condensate from the microscopic 
viewpoint. The Hamiltonian in a rotating frame 
with the angular frequency r2 rot is given by Ti = 

^(r). The creation and 

annihilation operators of the Bose particle are and 
which are decomposed into ^(r) = <f>(r) + V>(r). The 

single-particle Hamiltonian is given as h(r) = — ^2 

/i + V(r) — fi ro t • (r x p) with the confining potential 
V(r) and the chemical potential fi. The last term in Tt 
describes the interaction between the particles, which is 
classified into the resulting eight terms according to the 
noncondensate part ip, by using the decomposition of the 
field operator obtained above p6[ . 

The quadratic Hamiltonian may be diagonalized by the 
usual Bogoliubov transformation, 




This diagonahzation leads to the following conditions. 
First we impose the condition on the condensate wave 
function <p(r) as \h(r) + g\<fi(r)\ 2 ] <p(r) = 0, which is 
the so-called Gross-Pitaevskii (GP) equation. The 
Bogoliubov-de Gennes (BdG) equation for the quasipar- 
ticle is given in terms of the eigenfunctions w q and i> q 

{h(r) + 2 3 |0(r)| 2 } Wq (r) - # 2 (i> q (r) = £ q u q (r), 
{h(v) + 2<?|0(r)|> q (r) - ff 0* 2 (r)u q (r) = -e q » q (r), 1 



where q is the quantum number for the eigenstate. It 
is well known that the eigenvalues e q correspond to the 
normal modes of the condensate and the values are good 
agreement with experiment for systems with negligible 
thermal cloud (ljj. 

We consider the collective mode by the linear re- 
sponse theory fUSf , In the presence of an exter- 
nal field coupled to the density (^^r^ the linear re- 
sponse is characterized by a retarded correlation function 
D R {vv',t-t') = -i{[^{Tt)^{Tt),^{T't')^{r't')])9{t- 
t'), where (• • •} denotes the thermal average. In or- 
der to calculate D R directly, it is convenient to intro- 
duce a Matsubara correlation function X>(rr',r — r') = 
-([*t( rr )^( rr )^rt( r V)f (rV)]>. The Fourier coeffi- 
cient T>(rr' ,i£l n ) can be related to the retarded corre- 
lation function, D R (rr\ us), by using the analytic contin- 
uation id n — ► u> + ir\. Here r\ is a positive infinitesimal 
constant and we use rj = 0.005uv in our calculation. In 
the dilute Bose system, the Matsubara correlation func- 
tion is characterized by a matrix 

V(rv>, iSl n ) ~ - ( r (r), 0(r) ) Q(rr', ifi n ) ( ^ ) ,(3) 

where Q is a 2 x 2 matrix renormalized Green's function 
defined as £(rr,rV) = -(T r [i(rr)it(rV)]) with the 
matrix operator Jv = {ip' , ip) . Using the Beliaev-Dyson 
equation, we have Q = (Gq 1 — 0, |l^] with the bare 
Green's function Qq which is constructed from quasipar- 
ticles and the canonical transformation Eqs. ([!]) and (||). 

The self energy is expressed, by the second order per- 
turbation theory on g, as S = T,^ + & 2 \ The first order 
term & x \ which contributes to the energy shift, is 

E^WiA) = g J <M/ qi (r) ( «00 ^ J7 qi(r); (4) 

where the noncondensate density p(r) = 
E q [l u q( r )| 2 /( £ q) + l w q( r )| 2 {/( £ q) + 1}]) the anomalous 
average «(r) = - £) q u q( r )^ q ( r ){ 2 /( £ q) + 1}: and tne 
Bose function /(e q ) = [exp(/3e q ) — with (3 = l/fceT. 
The second order term X^ 2 ) is given by 
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which determines the excitation and its damping. In 
Eq. (||), the first term describes Landau processes: q[ + 
q' — > q — > qi + q', while the second and third terms 
correspond to Beliaev processes: q^ — ► q + q' — > qi. 
The matrix elements A and B are the overlap integrals 
of the condensate and the excitation amplitudes. The 
frequency spectrum of the collective modes is given as 



S(u) 



J dvj dv'F^{v)lmD R {vv'^)F Qe {v). 



(6) 



In the case of the surface mode (Q r = 0), the general 
excitation operator is given as Fq b = r^^e 1 ^" 9 

We take up the quadrupole excitation experiment on 
87 Rb atoms by Bretin et al. M . Assuming their long-cigar 
system as a cylinder, we introduce a cylindrical coordi- 
nate: r = (r,6,z). The quantum number q = (qg,q z ,q r ) 
may take the following values: qg = 0, ±1,±2, q z = 
0, ±2ir/L, ±47r/L, and q r = 0, 1, 2, where L is the 
period of the length along the z-axis and we take 
L = 15/im. The linear density along the z-axis can be es- 
timated as n z — 6 x 10 9 /m with the radial trap frequency 
w r /27r = 97.3Hz (T c ~ 200nK). At finite temperatures, 
following the procedure by Rusch et al. |M, we use the 
chemical potential /i to fix the total number. 

In Fig. 2(a) we show S(u>) in the frequency region near 
the quadrupole excitations u>±2, where we used the pa- 
rameters appropriate for the experiment by Bretin et al. 
M. The two resonance peaks are seen from it. The 
peak around u) — 1.55u; r corresponds to the main res- 
onance of the quadrupole mode qg = +2, while the 
other mode qg = —2 has a peak around ui = 1.3w r . 
These resonance frequencies correspond to the observa- 
tions (w +2 /27r = 159.5 ± 1.0Hz ~ 1.6uv/2tt, w_ 2 /2tt = 
116.8Hz ~ 1.2w. r /2vr). It is evident from Fig. 2(a) that 
the resonance width is wider in qg — —2 than that in 
qg = +2 by a factor two. 

Since the damping of the collective mode comes 
from the imaginary part of the self energy X^ 2 ) in 
Eq. (^]) and the dominant matrix element is E^, we 
discuss the damping constant defined as T qg (u>) — 

— ImE-j^ (q, q, Lj)\q r=qz =o. The damping constant for 
qg = —2 and qg — +2 are, respectively, shown in Fig. 2(b) 
and (c). For the quadrupole mode qg = +2, only the Lan- 
dau process (+2,0,0) + (qg,q z ,q r ) -> (qg + 2,q z ,q' r ) is 
active as in the vortex-free case. Thus the damping con- 
stant for qg = +2 is nearly same for both single-vortex 
and vortex-free cases. This is exactly seen by Bretin et al. 
||. On the other hand, the Beliaev process exclusively 
dominates the damping of the qg = — 2 mode because, as 
mentioned before, (—2, 0, 0) can decay into the two dipole 
modes w_i by conserving the angular momentum and the 
total energy. As is seen from Fig. 2(b), r_ 2 (w ~ 1.3co r ) 
for qg = — 2 at T = shows a single peak due to this 
particular Beliaev process while r +2 (w ~ 1.55ay) for 
qg = +2 shows no prominent peak (Fig. 2(c)). As T 
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FIG. 2: Spectrum of two quadrupole modes qe — ±2 at 
T — 15nK for rj — 0.005ij; r (a). The damping constant of the 
quadrupole mode qe = —2 (b) and qg — +2 (c) at T — 0, 5, 
10, and 15nK. The peak around uj — 1.3u> r corresponds to the 
Beliaev process (q e = -2, q z = 0) -> (-1, +2) + (-1, -2). 



increases, the Landau processes becomes effective, giving 
rise to the satellite peaks in addition to the main peak 
as shown in Fig. 2(b). We can see that at T ~ 15nK 
the ratio of two damping constant T±2 at each resonance 
is r_ 2 /r +2 ~ 2, which agrees with the observation M. 
However, this ratio depends sensitively on the tempera- 
ture; As T decreases, the Landau process quickly vanish, 
leaving only the particular Beliaev process active. This 
should be checked experimentally. 

The created modes via the Beliaev process of the 
counter rotating quadrupole mode qg — —2 can be iden- 
tified to the so-called anomalous mode |l2|, |l3) , or Kelvin 
mode characterized by the quantum number qg = — 1, 
q r = 0, and ±q z — Att/L as shown in Fig. 1(b). This 
dipole mode is counter-rotating to the vortex flow di- 
rection. Since this mode has the zero-relative angu- 
lar momentum to the condensate, the wave function 



u(r, z) 



= _i(r)e l<?zZ does not vanish at the center 



r = and localized at the core while all other wave 
functions with qg =/= — 1 vanish at r = 0. In Fig. 3 we 
display the condensate \<f)(r, z)\ and the wave function of 
the Kelvin mode with q z = dzAir/L, which is created by 
decaying of the counter-rotating quadrupole mode. One 
can see that this anomalous mode with a negative eigen- 



FIG. 3: Spatial profiles of the condensate \cj>(r, z)\ (a) and the 
Kelvin mode \u(r, z)\ = \u qe --\{r) cos(g z z)| with q z = ±4n/L 
(b). 



value at q z — 0, which is first identified theoretically 
Isoshima and Machida E3] and Dodd et al. is lo- 

calized within the core region along the radial direction 
whose scale is the coherence length. The characteristic 
wave number q z = Att/L approximately corresponds to 
the observation H. These particular features, the local- 
ization around the core in the radial direction and prop- 
agation along the vortex line direction, are exactly what 
Bretin et al. || have detected. Since there is no core 
localized mode other than this qg = — 1 mode [^2|, we 
conclude that a pair of these anomalous modes with q z 
and — q z are created. Thus the Kelvin wave, or the wave 
motion of the vortex line propagating along z-axis, is now 
identified and imaged by their experiment (see Fig. 3 in 
Ref. §). 

In the following we consider several situations to help 
identifying the Kelvin mode: Let us first consider the ex- 
ternal rotation effect on the decay process. Under the ex- 
ternal rotation frequency f2 ro t, ^±2(^rot) = ^>±i T 251 ro t 
as is seen from Eq. (0) ||. Accordingly, the spectral 
response function S(uj) shown in Fig. 4 exhibits; (1) 
The two resonances switch their positions at around 
Q ro t — 0.05w r , and (2) the resonance widths become com- 
parable as f2 ro t increases. It is due to the suppression of 
the Beliaev decay, because under the rotation the pop- 
ulation of the possible modes decreases according to the 
rule w 9fl (f2 r ot) = Lo qe (0) - qe^vot- 

The finite temperature affects both Beliaev and Lan- 
dau processes. As increasing T, w_i(T) is known to be 
larger while oj_2(T) is relatively independent of T ex- 
cept the region near T c Thus q z (T) becomes small 
with T and above some critical temperature the Beliaev 
process — > u>-i(q z ) + uj-i(—q z ) ceases to exist, un- 
abling to create the Kelvin modes as the decay channel of 
Generally the Landau process for both uj±2 modes 
becomes important as T increases because thermally ex- 
cited modes become more available. 

In Bretin et al. experiment [^| the counter-rotating 
mode is used to excite the Kelvin mode via the Be- 
liaev process. It is also possible to use the monopole 
mode luq to create the Kelvin mode, namely ujq — > 
u>+i(q z ) + uj-i{—q z ). This provides yet another spectro- 
scopic method to analyze the Kelvin mode. 
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FIG. 4: The external rotation effect on spectrum of the co- 
and counter-rotating quadrupole modes at 15nK. 



In conclusion, we have demonstrated that the Kelvin 
mode as a propagation wave along a vortex line can be ex- 
cited via Beliaev decay processes for the counter-rotating 
quadrupole mode. It enables us to understand the exper- 
iment by Bretin et al. j^] successfully, namely, their res- 
onance position and width, and to predict the external 
rotation and temperature effects. We have also shown 
that utilizing these decay channels provides a novel spec- 
troscopic method for low-lying Bosonic excitations in a 
vortex, in particular, the unexplored Kelvin mode. 

The authors thank P. Rosenbusch for communicating 
their results prior to publication. 



[1] 
[2] 

[3] 

[4] 
[5] 
[6] 
\<] 



[8] 

[9] 
[10] 

[11] 
[12] 

[13] 
[14] 



F. Dalfovo et al, Rev. Mod. Phys. 71, 463 (1997). 
A.L. Fetter and A. A. Svidzinsky, J. Phys.: Condens. 
Matter 13, R135 (2001). 

M.R. Matthews et al, Phys. Rev. Lett. 83, 2498 (1999). 
K.W. Madison et al, Phys. Rev. Lett, 84, 806 (2000). 
A.E. Leanhardt et al, Phys. Rev. Lett. 89, 190403 (2002). 
JR. Abo-Shaeer et al, Science 292, 476 (2001). 
In high T c cuprates the magnetic field induced vortex 
state is one of the focused themes. See for example, S. 
Sachdev and S.-C. Zhang, Science 295, 452 (2002). 
R.J. Donnelly, Quantized Vortices in Helium II (Cam- 
bridge Universit y Press, Cambridge , 1991). 



V. Bretin et al, cond-mat/0211101 



[15] 

[16] 

[17] 
[18] 



S.T. Beliaev, Sov. Phys. JETP 7, 289 (1958). 

P.O. Fedichev et al, Phys. Rev. Lett. 80, 2269 (1998). 

T. Isoshima and K. Machida, J. Phys. Soc. Jpn. 66, 3502 

(1997). 

R.J. Dodd et al, Phys. Rev. A 56, 587 (1997). 
It is generally accepted that it is negative with the 
positive norm conserved wave function. Virtanen et al. 
[S.M.M. Virtanen et al, Phys. Rev. Lett. 86, 2704 (2001)] 
claim that it could become almost zero even at zero- 
temperature. This point is still under discussion. See for 
review, Ref. j|. 

M. Rusch et al, Phys. Rev. Lett. 85, 4844 (2000); S.A. 
Morgan, J. Phys. B 33, 3825 (2000). 
A.L. Fetter and J.D. Walecka, Quantum Theory of Many- 
Particle Systems (McGraw-Hill, New York, 1971). 
S. Stringari, Phys. Rev. Lett. 77, 2360 (1996). 
T. Isoshima and K. Machida, Phys. Rev. A 59, 2203 (1999). 



